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Abstract: The abelian monopoles in lattice gluodynamics are associated with continuum 
monopoles blocked to the lattice. This association allows to predict the lattice monopole 
action and density of the (squared) monopole charges from a continuum monopole model. 
The method is applied to the static monopoles in high temperature gluodynamics. We 
show that the numerical data both for the density and the action of the lattice monopoles 
can be described in terms of a Coulomb gas of abelian monopoles in the continuum. 



o 

\SJ 

o 
o 

00 

o 
^1 



Keywords: [Nonperturbative Effects, Confinement, Lattice QCD 



© SISSA/ISAS 2008 



http: //jhep. sissa. it/ar chive/papers/ jhep092003027/jhep092003027 .pdf 



Contents 



[T]. Introduction |l| 

|2[ Lattice monopoles from continuum monopoles ^ 

|2.1| Lattice monopole action || 

|2.1.1| Action for large lattice monopoles ^ 

|2.1.2| Action for small lattice monopoles 10 

Validity of the gaussian approximation 10 

Lattice monopole density 12 



|3|. Numerical results 14 
|3.1| Details of simulations 14 



3.2 Lattice monopole action 17 



3.3 Lattice monopole density 18 



3.4 Check of the Coulomb gas picture for continuum monopoles 19 



EL Conclusions 20 



1. Introduction 

One of the most popular approaches to the problem of confinement of quarks in QCD is the 
so-called dual superconductor mechanism [Q]. The key role in this approach is played by 
abelian monopoles which are identified with the help of the abelian projection method Q. 
The basic idea behind the abelian projections is to fix partially the non-abelian gauge 
symmetry up to an abelian subgroup. For SU(iV) gauge theories the residual abelian 
symmetry group is compact since the original non-abelian group is compact as well. The 
abelian monopoles arise naturally due to the compactness of the residual gauge subgroup. 

In the dual superconductor mechanism the abelian monopoles are considered as effec- 
tive degrees of freedom which are responsible for confinement of quarks. According to the 
numerical results Q the monopoles are condensed in the low temperature (confinement) 
phase. The condensation of the monopoles leads to formation of the chromoelectric string 
due to the dual Meissner effect. As a result the fundamental sources of chromoelectric 
field, quarks, are confined by the string. The importance of the abelian monopoles is also 
stressed by the existence Q of the abelian dominance phenomena which was first observed 
in the lattice SU(2) gluodynamics: the monopoles in the so-called Maximal abelian projec- 
tion Q make a dominant contribution to the zero temperature string tension (for a review, 
see ref. ||). 



In the deconfinement phase (high temperatures) the monopoles are not condensed and 
the quarks are liberated. This does not mean, however, that monopoles do not play a 
role in non-perturbative physics. It is known that in the deconfinement phase the vacuum 
is dominated by static monopoles (which run along the "temperature" direction in the 
euclidean theory) while monopoles running in spatial directions are suppressed. The static 
monopoles should contribute to the "spatial string tension" — a coefficient in front of the 
area term of large spatial Wilson loops. And, according to numerical simulations in the 
deconfinement phase @ , the monopoles make a dominant contribution to the spatial string 
tension. 1 

In this paper we investigate the physics of the lattice monopole currents in the SU(2) 
gauge theory at high temperatures. The properties of the monopoles are effectively three 
dimensional due to static nature of the currents. However, the problem which is discussed 
in this article is quite general and it can not be ascribed only to the three dimensional case 
which is considered here as a useful and physically relevant example. 

The lattice monopoles in Monte Carlo simulations are detected by the standard De- 
Grand- Toussaint construction || which identifies the magnetic monopoles by measuring 
the magnetic flux coming out of lattice 3-D cells (cubes). The properties of thus-measured 
lattice monopoles should obviously depend on the physical size, b, of the lattice cell (below 
we call these lattice objects as "lattice monopoles of the size 6"). To get the continuum 
properties of the monopoles one should send the size of the lattice cells to zero, 6^0, what 
is usually a difficult numerical problem. Here we propose a method — which can be called 
as "blocking of monopoles from continuum to the lattice" — to identify the couplings of a 
continuum model for monopoles using the numerical results obtained on the lattice with a 
finite spacing. 

Each lattice 3D cell can be considered as a "detector" which measures the magnetic 
charges of the continuum monopoles travelling through the lattice: if the continuum mo- 
nopole is located inside a lattice 3-D cell then the DeGrand- Toussaint method detects a 
lattice monopole. From this point of view the limit b — > is needed to detect accurately 
the position of the continuum monopole with the help of the lattice detector. If the size of 
the lattice cell is finite then two or more continuum monopoles may be placed inside the 
cell. The fluctuations of the monopole charges of the lattice cells should be described by a 
continuum model. As a result, the lattice observables — such as the vacuum expectation 
value of the lattice monopole density — should carry information about dynamics of the 
continuum monopoles. The observables should depend not only on the size of the lattice 
cell, b, but also on features of the continuum model which describes the monopole dynam- 
ics. We take the high temperature gluodynamics as a simple example and show that the 
density and action of the static lattice monopoles of various sizes b can self-consistently be 
described by a continuum Coulomb gas model. 

Our approach resembles the idea of the blocking of the continuum fields to the lat- 



tice [ 10 1 — This method allows to construct perfect actions and operators in various 



1 Note that in another approach to the "spatial confinement" problem, the area law for the spatial Wilson 
loops is suggested to be caused by magnetic thermal quasi-particles discussed in ref. |B|. It seems plausible 
that the continuum abelian monopoles may be correlated with these magnetic excitations. 



field theories. In our paper we are blocking continuum monopoles what is ideologically 



similar to the blocking of a topological charge |ll|] and a fermionic current [10, 12] to the 
lattice. The blocking of the fields, however, leads to a non-integer lattice magnetic current 
contrary to the quantized lattice definition of the monopole charge ||]. Thus the blocking 
of topological defects seems to be more suitable for investigation of the lattice monopole 
charges. There is also a similarity of our approach with the blocking of the lattice monopole 
degrees of freedom from fine to coarser lattices |l4j]. This method allows to define a perfect 
monopole action independent on the spacing of the fine lattice. 

To avoid misunderstanding we would like to stress from the very beginning the differ- 
ence between various lattice monopole sizes b. As we mentioned above, we call the size 
of the lattice 3D cell — used for detection of the monopoles — as "the size of the lattice 
monopole". This should be distinguished from the physical radius, ro, of the monopole 
core [15] which is obviously independent on the size of the lattice "detector". In this paper 



we disregard the existence of the monopole core and consider the continuum monopoles as 
point-like objects. Thus, we get the finite-sized lattice monopoles by blocking the point-like 
monopoles from continuum to the lattice. 

There is another motivation for our study. Various properties of the lattice monopoles 
such as the action and density of the monopole currents, depend on the physical size of 



the lattice monopole [16, [17]]. This effect is observed both at zero 

M an d finite HI tem- 
peratures. Usually the monopole density is calculated numerically for so called elementary 
monopoles the size of which coincide with the lattice spacing a. Then the monopole quan- 
tities in continuum limit (e.g., continuum monopole density) are usually associated with 
the limit a — > 0. However, since the monopole in the gluodynamics has a core of the finite 
size [|l5|] the continuum monopole may simply "drop" through the lattice cells if a <C ro. 
Moreover, the scale of the ultraviolet lattice artifacts coincides with the lattice UV cut-off 
a, which, in turn, equal to the size of the elementary monopoles. Therefore the lattice 
elementary may be affected by the lattice UV artifacts. Both these reasons imply that (at 
least, naively) there is a possibility that the continuum monopole density may not cor- 
rectly be calculated as the elementary lattice monopole density in the limit a — > 0. One of 
advantages of our approach is that it allows to relate the lattice monopole quantities to the 
corresponding physical quantities in continuum limit at the finite lattice spacing making 
irrelevant the problems noticed above. 

The structure of the paper is as follows. In section [2] we derive the lattice monopole 
action and the density of squared magnetic charges. The basic assumption behind the 
derivation is that the dynamics of the continuum monopoles — which are blocked to the 
lattice — is described by the Coulomb gas model. In section || we present the numerical 
results both for the monopole action and for the density of squared magnetic charges. We 
show that these quantities are in a good agreement with the predictions of section ^. The 
comparison of the analytical and numerical results allows us to calculate the product of 
the abelian magnetic screening mass and the monopole density in the continuum model. 
At the end of section ||| we check the self-consistency of our results as well as validity of the 
Coulomb gas description for the static monopole currents. Our conclusion is presented in 
the last section El 



2. Lattice monopoles from continuum monopoles 



In this section we consider the blocking of the continuum monopoles to the lattice in three 
space-time dimensions. Let us consider a lattice with a finite lattice spacing b which is 
embedded in the continuum space-time. The cells of the lattice are defined as follows: 
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b[ Si - < Xi < b(^Si + - 



1,2,3 



(2.1) 



where s, is the lattice dimensionless coordinate and Xi corresponds to the continuum co- 
ordinate. 

Each lattice cell, C s , detects the total magnetic charge, k s , of "continuum" monopoles 
inside it: 



d 3 x p(x) 



p(x) = Y,<l a 5 ia) (3 



.(a)> 



(2.2) 



c s 



where p{x) is the density of the continuum monopoles, x a and q a is the position and the 
charge (in units of a fundamental magnetic charge, qm) of a th continuum monopole. In 
three dimensions the monopoles are instanton-like objects and the monopole trajectories 
have zero dimensionality (points). It is worth stressing the difference between continuum 
and lattice monopoles: the continuum monopoles are fundamental point-like objects in 
this approach while the lattice monopoles correspond to the lattice cells with non-zero 
total magnetic charges of continuum monopoles located inside appropriate cells. 

According to definitions ([2.2;), the lattice monopole charge shares similar properties to 
the continuum monopole charge. The lattice monopole charge k s is quantized, k s £ Z, and 
conserved in the three-dimensional sense: 



sGA 



d 3 x p{x) 



0. 



(2.3) 



if the continuum charge is conserved. Here A and V denote the lattice and continuum 
volume occupied by the lattice, respectively. In other words, the total magnetic charge of 
the lattice monopole configuration is zero on a finite lattice with periodic boundary which 
is considered in this paper. 

Suppose that the dynamics of the continuum monopoles is governed by the 3-D Cou- 
lomb gas model: 
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(2.4) 



a, 6=1 



The Coulomb interaction in eq. ( p.4| ) is represented by the inverse laplacian D: —dfD(x) = 
8^(x). The density of the continuum monopole is assumed to be low. The continuum 
monopole charges therefore are restricted by the condition \q a \ < 1 which means that these 
monopoles do not overlap. The average continuum monopole density p is controlled by the 
fugacity parameter £. 
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Note that the model ( |2.4D does not exist without properly defined ultraviolet cut-off. 
Indeed, the self-energy of the point-like monopoles is a linearly divergent function due to 
the infinite nature of the Coulomb interaction. This divergence — which is not explicitly 
present in eq. ( |2.4[ ) — should appear manifestly in a field representation of the Coulomb gas 
model which we need below (see, for example, the sin-Gordon representation for the lattice 
monopole action, eq. ( |2.10D ). The renormalized fugacity is, Cren = C " ex -P{9m/(& 7r r o)}j 
where the ro is the ultraviolet cut-off. In our case this cut-off is given by the size of the 



monopole core which is of the order of 0.05 fm at zero temperature [15|. For simplicity we 
omit the subscript "ren" in the renormalized fugacity below. 

The magnetic charges in the Coulomb gas (^J|) are screened: at large distances the two- 
point charge correlation function is exponentially suppressed, (p(x)p(y)) ~ exp{— \x — y\/ 
Xd}- Here Xd is the Debye screening length 111 



X D 



9Myfp 



(2.5) 



Note that the three dimensional Debye screening length corresponds to a magnetic screen- 
ing in four dimensions. The density of the continuum monopoles in the leading order is 
related to fugacity as |l^] p = 2Q. Below we choose the vacuum expectation value of the 
continuum monopole density, p, and the Debye screening length, Xd, as suitable parameters 
of the continuum model (instead of qm and £) . 

We are interested in two basic quantities characterizing the lattice monopoles: the 
lattice monopole action S mon (k) and the v.e.v. of the squared magnetic charge, (k^). We 
study the quantity (kg) instead of the density, (|/c s |), since the analytical treatment of the 
density is difficult while both these quantities are equivalent philosophically. 

2.1 Lattice monopole action 

To derive the lattice monopole action we substitute the unity, 



i= e ( k >- f d 3 ^)) > 

fe(A) 6 ZseA \ c s J 



(2.6) 



into eq. fl2.4p . Here J]fc(A)ez = YlseA^2k s ez an< ^ stands here for the Kronecker symbol 
(i.e., lattice ^-function). We get 



Z 



~/\n 



fc(A)eziv=o 



x exp 



f\ J d 3 xW [v A h[v X x 

a=l q a =±lJ _ n 

— (diX(x)) 2 + ¥i x ) [ X{x) - d s(x)h s I k sh s \ , (2 

g M \ seA J seA J J 



where we have introduced two additional integrations over the continuum field % an d the 
compact lattice field h to represent the inverse laplacian in eq. (|2.4| ) and the Kronecker 
symbol in eq. ( |2.6D , respectively. The subscript A in T>\h indicates that the integration is 
going over the lattice fields h. The representative function of the s th lattice cell is denoted 
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as 0, 



1 , x G C s , 
, otherwise . 



Summing over the continuum monopole positions according to ref. [18], we get 

k(A)eZ 

where the lattice monopole action, S mon , is defined as follows: 



(2. 



(2.9) 



Jmon 



Dx / 2?A^exp< 



d 3 x 



1 



2& 



(^x(x)f - 



(2.10) 



2C cos x^) - d s(x)h s ) k s h s 



seA 



The functional integral over the field x i n the l.h.s. of this equation resembles the sin- 
Gordon model coupled to an external source. 

Note that the left hand side of eq. ( 2.1C| ) is invariant both under the global continuum 
transformations of the field x> an d under the lattice local transformations of the field h, 



x(x) — > x( x ) + 2vrn , n G Z 
h s — ► h s + 2irm s , m s G Z , 



(2.11) 



respectively. Due to the last invariance we extend below the integration over the lattice 
field h to infinite limits (this leads to appearance of an inessential factor in front of the 
partition function). 

Let us consider the lattice monopole action in the tree (or, gaussian) approximation. 
The validity of this approximation is discussed in section |2.2j . We expand the cosine function 
over the small fluctuations in x and h and we get in a leading order: 



^mon CO 



T>x / T>\hexp< 



d 3 x 



1 



2<4 



(diX{x)) + 



(2.12) 



V sGA 

The integration over field x leads to the following expression: 



where 



/ V A hex.p< -(b 3 ^2 h s F s> lih s ' + i ^ k s h s > , 
I s,s' seA ) 



Gs,s> = ^ j d 3 x j d 3 y6 s (x)D mD (x -y)e s ,(y) , 



(2.13) 

(2.14) 
(2.15) 



where D mD is the scalar propagator for a massive particle, (— df + m 2 )D m (x) = 5^(x), 
with the Debye mass m = mo = A^ 1 . Note that the lattice operators T and Q are 
dimensionless quantities. 

In eq. ( |2,13| ) the integration over the lattice field h can be performed straightforwardly 
giving us the action for the lattice monopoles: 



(2.16) 



S,S' 



Let us calculate the operator T on the infinite lattice. We represent the propagator 
D mD as an integral over continuum momentum and integrating over x and y coordinates 
in eq. (2.15): 



1 r d 3 P 1 

y s . s ' — TIT 



6 5 J (2ir) 3 p* + m 2 D 



n 

i=l 



2sin(6pi/2) 
Pi 



J(s-s',p)b 



(2.17) 



where (a, b) denotes the scalar product of the vector quantities. Changing in eq. ( pT7D the 
integration variable, p = q/b, and introducing the dimensionless mass, 



we get 

G s ,s' 



d 3 q 



i(s — s',q) 



fj. = b/X D , 



H{q,p) = 2 - 2 



n 

i=i 



sin(g i /2) 
ft/2 



(2.18) 



(2.19) 



Now we have to find the operator .F which is inverse to eq. ( 2.14 ). We represent F as 
follows: 

T K j = [ n d 3 uT(u) e i (»-»'.«) , (2.20) 



Substituting this equation and eq. ( |2.14| ) in the relation, 



s "eA 



we get the equation, 

d 3 u \ F(u) 



1 - /i 2 ^ H(u + 27rr, p) 

rGA 

where we have used the Poisson summation formula, 



J(s-s',u) _ Q ^ 



e *('">«-9) = (2vr) 3 ^5 (3) (g -u- 2vrr) , 
s"GA reA 



(2.21) 



(2.22) 



and then integrated over the momentum q. 
The solution to eq. ( 2.21| ) is 

1 



T{u) 



(27r)3l-/i 2 E r6 A^(« + 2 ^^) 



(2.23) 
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Substituting eq. ( [2.23 ) in eq. ( EDcD and using definition ( |2.19| ), we finally get 
d 3 u \ ^ ^ 4sin 2 (u i /2) A / 2 sin(uj-/2) \ 2 ] 



(2k)* 



f,t 



Lr(A)gZ t 

where we have used the relation 

E 



^ (u + 2k?) 2 + fi 2 *-^yuj + 



n 



2kt~ 



3= 



■ e 



i(s—s',u) 



(2.24) 



(u i + 2Kr i ) 2 [4 sin 2 (ui/2)] ' 



(2.25) 



The operator !F S s > coincides with a three-dimensional perfect propagator for a free 
scalar field space-time discussed in details in ref. [p4|]. 

The summation in eq. ( |2.24j ) over one of the integers ri can explicitly be done with the 
help of the relation, 

sinh A j (u + 2k? , fi) 



E 



1 



1 



where 



(u + 2k?) 2 + fi 2 2Xi(u + 2k?, fi) cosh Aj(u + 2k?, fi) — cos Ui 

1/2 



(2.26) 



E«J + 



3=1 



(2.27) 



Substituting eq. fljTgg ) into eq. (|2.24| ) we get the final expression for operator JF: 
d 3 u 



(27T)' 



e i(s-s',w) x 



(2.28) 



EE 

1=1 r 3 -6Z 
37^ 



2 sin 2 (uj /2) sinhAi(u + 27rf,^) -,4 4sin 2 (u A ./2) 

\i(u + 2k?, fi) cosh Xi(u + 2k?, fi) — cos U{ ^ + 2Kr^) 2 



The leading term in the lattice monopole action is defined by eqs.( 2.16| ), ( [2.280 . 

The finite-volume expression for the lattice monopole action can be obtained from 
eq. ( p. 28 ) by the standard substitution: 



11; 



2Kki 
~L~ 



k = 0,1,... Li - 1 



2k U 



Li-l 



1 ~* " 

rE- 



(2.29) 



ki=0 



where Li is the lattice size in z th direction. 

In the infinite-volume case the lattice operator T s s i depends only on the dimensionless 
quantity fi, eq. ( p. 18 ), which is the ratio of the lattice monopole size b and the Debye 
screening length, eq. (|2.5| ). As we will see below the form of the operator T is qualitatively 
different in the limits of small and large fi. Thus the Debye length Ad sets a scale for the 
lattice monopole size (or, better to say, for the size of the lattice cell) which characterizes 
different forms of the lattice monopole action. 
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2.1.1 Action for large lattice monopoles 

Let us consider the case of the large lattice monopoles, b 3> Xd, or, equivalently, /j>1. In 
this case eq. ( |2.28 ) can be simplified since functions Aj (eq. (|2.27 )) are large and sinhAj ~ 
coshAj ~ e Xt /2. Up to 0(e _/i ) corrections the operator ,F S s ' is given by 



d 3 u 



2sin 2 (uj/2) -j^j- 4sin 2 (u fc /2) 



i=l r^SZ 



■ e 



i(s—s',u) 



(2.30) 



In the case of large p, the function Aj in eq. (|2.27 ) is close to p for all r-j <C p/ (2ir). The 
deviation of A from p becomes substantial only for terms with large summation variable 
rj, Ti 3> p/(2tt). However such large r^'s are in any case suppressed as ~ Wj^i (2 7rr j)~ 2 
according to eq. (|2.30| ). Therefore the approximation A, = p would lead only to 0(p~ 2 ) 
corrections in eq. ( |2.30|) . Up to these corrections eq. ( 2.30 ) reads as follows: 



^, s '(A*»l) = f 



d 3 u 
(27p 



E- 2 (f)iI£ 



i=l 



4sin 2 (uj/2) 
(uj + 2'Krj) 2 



. p i{s-s',u) 



(2.31) 



In turn, eq. ( |2.3l| ) can be simplified with the help of the relation ( 2.25| ), giving: 

F S)S ,(p.^l) = 2pD S)S ,, (2.32) 
where D s s i is the lattice laplacian: 



d^u 
(2^)3 



4 E 



8111 



11 j 



J(s-s',u) 



(2.33) 



Thus for large sizes of the lattice cells, 6 > Ao, the lattice monopole action ( 2.16| ) is 
of the Coulomb-type: 



STol(k) = C(b) k * D *,s> K> + o( (^j V b » A D . 



(2.34) 



The pre-Coulomb term C can be expressed through the continuum density of the mono- 
poles, p, the Debye screening length A^> and the lattice spacing b as 



C(b) 



1 



1 



(2.35) 



XdP b 2 ' 

Here we used eq. ( |2.18 ). 

Thus in the limit of large lattice monopole sizes, b, the lattice monopole action is of the 
Coulomb form. The long-range nature of the Coulomb interaction in the effective lattice 
action and the screening effect in the underlying high temperature theory do not contradict 
each other as one may think. The analog of this situation in continuum can be represented 
by a similar Coulomb gas model ( |2.4[ ) which contain the long-range interaction part in the 



lagrangian. However, the interactions between monopoles lead to the screening effect in 
various correlators in this model. Obviously, the lattice Coulomb gas model ( [2.34 ) must 



also possess the screening effect despite the long range nature of the Coulomb term in the 
action. 

The coefficient in front of the Coulomb term, C(b), scales as C(b) ~ b~ 2 . Note that 
this form of scaling is a non-perturbative effect. Indeed, naively one could expect that this 
coefficient has to be proportional to the (squared) renormalized magnetic charge in three 
dimensions, C(b) ~ g 2 M . The magnetic charge is inversely proportional to the 3D electric 
charge, gu ~ 9e\d- ^ n ^ ne lading (tree) order of dimensional reduction formalism the 
electric charge is given by <?e,3d = <74d(A, T) ■ vT, where T is the temperature, g^D is the 
running AD charge of the SU(2) gluodynamics and A is the renormalization scale. As a 
result, C(b) ~ [gj D (A,T)T]~ . The size of the lattice monopole, b, may enter the above 
expression only in the form of the renormalization scale, A ~ 1/6. However, this would 
lead only to logarithmic 6-dependence of the coefficient C(b). Thus the b~ 2 dependence of 
the pre-Coulomb coefficient is clearly of a non-perturbative nature. 

2.1.2 Action for small lattice monopoles 

In the case of small lattice monopoles the leading term of the operator T for b <C Ad follows 
immediately from eq. ( [2.14 ): F s s i = 5 S s /. Thus the monopole action ( |2.16| ) becomes: 

S^tKk) = M(b) J2k 2 s + o(^-^j^j , 6 « Ad , (2.36) 
where the coefficient 

M(b) = 1 ■ i (2.37) 

plays a role of the lattice monopole mass. Indeed, the leading term in the action for N 
small lattice monopoles with the charges q = ±1 is S = M(b) ■ N. 

In summary, we have established that the monopole action depends on the ratio of the 
lattice cell and the continuum Debye screening mass, b/\o'- 



1 1 'v^ ( 2 - 38 ) 

— — • — • ) k s D s s , k s i H , b > A D , 

p\ D ¥ ^ 



where D s s < is the inverse laplacian on the lattice. Thus the leading contribution to the 
lattice monopole action is given by the mass (Coulomb) terms for small (large) lattice 
monopoles. 

2.2 Validity of the gaussian approximation 

As we will see below, from the point of view of numerical computations the most inter- 
esting and reliable case corresponds to the large-6 monopoles. Here we briefly discuss the 
validity of the gaussian approximation in this limit. The detailed derivation will be given 
elsewhere [19 |. 



We have truncated the cosine function in eq. ( |2.10| ) and performed the gaussian integra- 
tion over the continuum field x m ec l- Q2.12| ) to get the effective action in the leading order 
for the lattice field h s . This action appears under the exponential function in eq. ( 2.13; ). 
The integration over the lattice fields h led us to the Coulomb action ( 2.34j ) for the lattice 
monopoles in the large-6 limit. Let estimate the Coulomb action Q2.34 ) coming from the 



next order of the truncation of the cosine function in eq. (2.10). 

Expanding the cosine function in eq. ( |2.10| ) up to the fourth order and treating the 
correction as a perturbation we arrive to the following correction to the monopole action, 




9 a (x)h 




(2.39) 



where we used the superscript "(4)" because, as we show below, this action contains, 0(k A ) 
monopole-monopole interactions. The quantum average ( 2.39| ) is taken in the partition 



function ( 2.12 ) containing the monopole current k as a source. 

Integrating the fields x an d h, and neglecting independent of k terms we arrive to the 
following A; 4 -correction to the monopole action: 
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where 
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(2.41) 



The operator Y is a converging function both in the limits of large and small q. 

According to eq. ( 2.32j ) the operator T is proportional to the first power of b in the 
limit n = m£)b S> 1. The estimation of the operator Y, eq. ( 2.41|) , is difficult even in 
the limit of large b. However, one can notice that Y contains the dependence on b in the 
denominators under the integrals over q, and, similarly to the operator J- -1 , the operator 
Y should not be a rising function of the scale b. Estimating Y as being of the order of 



unity, we get: |S'mon(^)| ^ g\.j /(48C6 5 )(3(fc 4 ). Thus, in the limit b — > oo the action 5mon 



:(4) 



provides a small correction to the Coulomb action ( p. 34 ), which is proportional to b 2 . 
These results are confirmed by the fact that in the limit of large monopoles k 4 correction 



to the monopole action is small compared to the quadratic terms, ref. [17|. Moreover, as we 
will see below, our numerical calculations shows that the quadratic part of the monopole 
action can be described by the Coulomb interactions with a high accuracy. 



2.3 Lattice monopole density 

In this section we discuss the dependence of the density of the extended lattice monopoles on 
the monopole size b. The simplest quantity characterizing the monopoles is the monopole 
density pi a tt(b) measured in the lattice units: 



Piatt (6) = -p (J2\k 



(2.42) 



\seA 



where L is the lattice size in units of b. However, the treatment of eq. ( [2.42 ) from the 
analytical point of view is more difficult than that of the squared density of the squared 
lattice monopole charges: 



(k 2 (b)) = j- 3 




(2.43) 



Below we discuss the quantity ( 2.43| ) which has a similar physical meaning to the lattice 
monopole density fl2.42|) . 

Both in eq. (|2.42j ) and eq. ( 2.43 ) the value of k s is equal to the total magnetic charge 
of the continuum monopoles which are placed inside the corresponding cube of the volume 
b 3 . Obviously the lattice density ( 2.42] ) does not give the correct density of the continuum 
monopoles in general since if the lattice monopole size is large enough then the oppositely 
charged continuum monopoles cancel fields of each other inside the cell C s . Moreover, the 
dependence of the densities ( |2,42| ) and ( 2.43| ) on the lattice monopole size b must reflect the 
dynamics of the continuum monopoles. Indeed, one can expect that the functions piatt(^) 
for the Coulomb monopole gas and for the random monopole ensembles should differ from 
each other since the inter-monopole correlations are absent in the latter case contrary to 
the former one. As we will see below the situation is similar to the lattice monopole action 
discussed in the previous section. 

Using eq. (U) the lattice density can be written in the continuum theory as 

follows: 

(k 2 }= I d 3 x / d 3 y(p(x)p(y)}, (2.44) 

J Cs J Cs 

where the lattice site s is fixed and the average is taken in the Coulomb gas of the magnetic 
monopoles described by the partition function ( |2.4| ). We assume the validity of the dilute 
gas approximation for the continuum monopoles. 

The correlator of the continuum monopole densities, {p(x)p(y)), is well known from 



ref. [ 18 1 . Introducing the source for the magnetic monopole density, J, eq. ( 2.44 ) can be 
rewritten as follows: 



(k 2 ) 



Cs dX J Cs dy SJ(x)5J(y) V 



( exp <^ i / d 6 zp{z) J{z) 



(2.45) 



j=o 



Then we repeat the transformations in the previous section which led us to eq. ( 2.12j ). 
Integrating over quadratic fluctuations of the field x we get in the leading order 
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(p(x) p{y)) = p <5 (3) (x-y)- m 2 D D mD (x - y) 



(2.46) 
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Substituting eq. Q2.46 ) in eq. (|2,44j ) and taking the integrals over the cell C s we get 

(k 2 ) = pb 3 P(p), (2.47) 



where p is the dimensionless Debye mass, eq. ( 2.18] ) and in the thermodynamical limit the 
function P is 

d 3 q 



P(ja) = (F-\ (ji) = l-^[ 7^ H(q, p) . 

J (2vr) 



(2.48) 



Here the inverse matrix J- S%s i is given by eq. fl2.14|) and the function H is defined in eq. ( 2,19j ). 
The finite- volume analog of eq. ( |2.48| ) can be easily written using the substitution ( 2.29j ), 
Note that eq. ( \2.4 r ' ), (2.4i ) establish a direct relation between the density of the squared 



lattice monopole charges and the lattice monopole action, eq. (2.16). 

It is interesting to study the scaling laws of the lattice monopole densities for large 
and small monopoles. In the limiting cases the behaviour of the function P is as follows: 

( Ci fx' 1 + 0(fi- 3 ) , p 
1 - C 2 p 2 + 0(p 4 ) , p 



where the constants Ci ; 2 are 

Ci » 2.94 . 



CXD , 

0, 



(2.49) 



Co 



d 3 q 

W 3 



H(q,0) « 0.148. 



(2.50) 



Substituting the asymptotic functions ( |2.49| ) in eq. ( 2.47| ) we get the scaling laws (in 
physical units) for the density of the squared lattice monopole charge: 



(k 2 



C lP X D b 2 
pb 3 ■ 



i + o (x D /by 



i + c 2P (b/x D ) 2 +o({b/x D ) 4 



(2.51) 



b-> A D , 
b < Ad • 

One can note some interesting properties of the density of the squared magnetic charge, 



eq. (HI). 

(i) The dependence of the lattice monopole charge squared on the monopole size b is 
always polynomial, b n . In particular, this property is remarkable in the large-6 region: 
one may expect that the screening of the monopole charges in the magnetic monopole 
plasma would lead to an exponential behaviour, e - cons t- m D fe ) which is not the case. 

(ii) As in the case of the lattice monopole action, the power n of the leading scaling law 
depends crucially on the value of the ratio of the lattice monopole size and the Debye 
screening length. 

(hi) The proportionality of the density (k 2 ) to b 2 in large-6 region has a simple explana- 
tion. In a random gas of continuum monopoles we would get (k 2 ) ~ pb 3 since the 
monopoles are not correlated with each other. In the Coulomb gas the monopoles 
are correlated and moreover, screened. Therefore the monopoles separated from the 
boundary of the cell by the distance larger than Ad ; do not contribute to (A; 2 ). Con- 
sequently, the b 3 proportionality for the random gas turns into Xpb 2 in the Coulomb 
gas and we get (k 2 ) ~ pXj^b 2 . The coefficient of proportionality is of a geometrical 
origin. 
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(iv) In the small b region the density of the squared lattice monopole charges is equal to 
the density of the continuum monopoles times the volume of the cell. This is natural, 
since the smaller volume of the lattice cell, 6 3 , the smaller chance for two continuum 
monopoles to be located at the same cell. Therefore each cell predominantly contains 
not more that one continuum monopole, which leads to the relation kg = \k s \ = 0, 1. 
As a result we get (k 2 ) — ► pi a tt(&) —> pb 3 in the limit b — > 0. 

Closing this section we mention interesting relations between the density of the small- 
and large- sized lattice monopoles and the coefficients in front of, respectively, the mass 
and the Coulomb terms of the monopole action, eqs,( |2.38 ), ( 2.51] ) : 



c{b) = wk)' 6>>Ad ' and M{b) = AWW)' b<<x °- (2 - 52) 

All these results are obtained in the gaussian approximation. The results on next to 
the leading order correction to the action and the monopole density will be published 



elsewhere [19]. 



3. Numerical results 
3.1 Details of simulations 

In order to get perfect lattice monopole action and density we perform numerically block- 
spin transformations for the lattice monopole currents. The original model is define on the 
fine lattice with the lattice spacing a and after the blockspin transformation, the renor- 
malized lattice spacing becomes b = na, where n is the number of steps of blockspin 
transformations. The continuum limit is taken as the limit a — > and n — ► oo for a fixed 
physical scale b. 

Finite temperature system possesses a periodic boundary condition for time direction 
and the physical length in the time direction is limited to less than 1/T. In this case it is 
useful to introduce anisotropic lattices. In the space direction, we perform the blockspin 
transformation and the continuum limit is taken as a s — > and n s — > oo for a fixed physical 
scale b = n s a s . Here a s is the lattice spacing in the space directions and n s is the blockspin 
factor. In the time direction, the continuum limit is taken as at — > and Nf — ► oo for a 
fixed temperature T = l/{Ntdt). Here at is the lattice spacing in the time direction and 
Nt is the number of lattice sites for the time direction. After taking the continuum limit, 
we finally get the effective lattice monopole action which depends on the physical scale b 
and the temperature T. 

In this paper the numerical procedure to generate the field configurations is identical 
to the one used ref. [17]. While referring an interested reader to the above paper, here 
we mention for completeness briefly the basic points of our numerical procedure. The 
anisotropic Wilson action for pure four-dimensional SU(2) QCD is written as 

P^(s) = \ Tr [l ~ U^ S )U„(s + £)£/t( s + u)ul(s)] + h.c. (3.1) 



p 


7 


a s 


a t 


P 


7 


a s 


at 


2.470 


2.841 


0.250 


0.075 


2.565 


2.152 


0.180 


0.075 


2.500 


2.615 


0.225 


0.075 


2.573 


2.098 


0.175 


0.075 


2.533 


2.354 


0.200 


0.075 


2.581 


2.042 


0.170 


0.075 


2.548 


2.256 


0.190 


0.075 


2.598 


1.927 


0.160 


0.075 



Table 1: Parameters /3, 7 and corresponding lattice spacings a s , a t in units of the zero-temperature 
string tension, y/<TT=o- 



If 7 = 1, the lattice is isotropic (a s = at). We can consider various lattice spacings a s 
and at by varying the parameters (3 and 7 in the action. The procedure to determine the 
relation between the lattice spacings a s , a t and the parameters f3, 7 is described in details 
in ref. fli7| ]. Using the parameters which are determined in ref. [17|, we simulate the pure 
four-dimensional SU(2) QCD on the lattice 48 3 x Lt with Lt = 12, 10, 8, 6, 4 corresponding 
to temperature T/T c = 1.6, 1.92, 2.4, 3.2, 4.8. The parameters used here are summarized in 
Table |]. The scaling behaviors for time-like lattice monopole action at above mentioned 
lattices show in large b region (T^ |. 

To study the abelian monopole dynamics we generate the thermalized non-abelian link 
fields {Ufj,(s)} and we perform abelian projection in the Maximally abelian (MA) gauge || 
for each SU(2) configuration. The MA gauge fixing condition is the maximization of the 
quantity R, 

maxR[U^} , R[U] = Tr £ [c^(a)<7 3 C^(s + £)<7 3 ] . (3.2) 

under the SU(2) gauge transformations, U — ► U^ 1 ' = The gauge fixing condition 

d3.2j ) is invariant under an abelian subgroup of the group of the SU(2) gauge transforma- 
tions. Thus the condition ( |3.2[ ) corresponds to the partial gauge fixing, SU(2) — ► U(l). 

After the MA gauge fixing, the abelian, {u^s)}, and non-abelian {U^^s)} link fields 
are separated, U^(s) = C^(s)u )J ,(s) , where 

C (s) -(V^WW A u(s) -( eldAS) ° ) (3 3) 

^ () "V c,{s) Vl-I^WlV' M() "V e-^J- (3 ' 3) 

The vector fields C^s) and u^(s) transform like a charged matter and, respectively, a 
gauge field under the residual U(l) symmetry. Next we define a lattice monopole current 
(DeGrand-Toussaint monopole) ||. Abelian plaquette variables 6^ v (s) are written as 

M s ) = M) + M + £)- M + i>) - 6 u (s) , (-4tt < M s ) < 47 • I 3 - 4 ) 
It is decomposed into two terms using integer variables n^ v (s): 

M s ) = M s ) + 2vrn^(s) , (-vr < M s ) ^ *")■ ( 3 - 5 ) 

Here 0^ u (s) is interpreted as an electromagnetic flux through the plaquette and n^ u (s) 
corresponds to the number of Dirac string piercing the plaquette. The lattice monopole 
current is defined as . 

M s ) = 2 e ^u P ad u np a (s + p) . (3.6) 

It satisfies the conservation law d'^k^(s) = 0. 



To study numerically the lattice monopoles of various lattice sizes we use the so-called 



extended monopole construction [20]. At zero temperature the extended monopoles can 
be defined in the symmetric way. They have the physical size b 3 where b = na. The 
charge of the n-blocked monopole is equal to the sum of the charges of the elementary 
lattice monopoles inside the n 3 lattice cell. At finite temperature the blocking of spatial 
and temporal currents should be done separately |L7]]: 

n s — 1 rit — 1 

K^a(s s , s 4 ) = ^2 ^2 ^ 4 ( HsSs + ( Us ~ ^ + ™ + I?' ntSA + ^ ' 

ij=0 1=0 
n a — l 

K4(s s ,Si) = ki (n s s s + ifi + ji> + //S,ra t s 4 + (n t - 1)) , (3.7) 
i,j,i=a 

where n s (nt) is the number of blocking steps in space (time) direction. 

We consider only the nt = 1 case since we are interested in high temperatures for which 
the monopoles are almost static. The lattice blocking is performed only in spatial directions, 
n s = 1 . . . 8, and we study only the static components K± among the 4Z) monopole currents 

(below we denote K4 as k.). At high temperature we disregard the spatial currents Ki 
since they are not interesting from the point of view of the long-range non-perturbative 
spatial physics. The size of the lattice monopoles is measured (unless otherwise specified) 
in terms of the zero temperature string tension, ctt=o- 

In order to figure out whether the continuum monopole currents can be described 
by the Coulomb gas model we should compare Monte Carlo results with the appropriate 
analytical predictions derived in the previous section. In principle this should allow us to 
obtain all independent parameters of the Coulomb gas (the continuum monopole density 
and the Debye screening length). However, the lattice monopoles of small sizes are largely 
affected by the lattice artifacts since in this case the number of blocking steps is small 
and possible magnetic charges of such monopoles are restricted due to peculiarities of 
the DeGrand-Toussaint definition ||. Moreover, the short-range interaction between the 
continuum monopoles should deviate from the simple Coulomb law due to a non-zero finite 
radius of the abelian monopole |jl5[ |. Thus in order to get reliable results we perform the 
comparison of the numerical data with analytical predictions for sufficiently large blocking 
steps only. The restriction to large-6 region allows us to calculate the product of the 
continuum monopole density and the Debye screening length, 

R(T) = X D (T) p(T) , (3.8) 

while a separate calculation of these quantities is not possible. Nevertheless the knowledge 
of quantity is enough to make a conclusion about realization of the Coulomb gas picture 



for static continuum monopole lines. In subsections |3.2| and 3.3 we get the quantity R 
from the lattice monopole action and density, respectively. To avoid confusion we denote 
R in these cases as R ac t and R p , correspondingly. In subsection 'iA we check the consis- 
tency of the obtained values of R ac t an d R p with each other and with the Coulomb gas 
picture. 



3.2 Lattice monopole action 



The lattice monopole action for the static mo- 
nopole currents, k s , at high temperatures was found 
numerically in ref. [|17]] using an inverse Monte- Carlo 
procedure. It turns out that the self-interaction of 
the temporal lattice currents can be successfully de- 
scribed by the quadratic monopole action: 
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(3.9) 
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Figure 1: The coupling constants, fi, of 
the lattice monopole action vs. the dis- 
tance, r at T — 1.6T C for various spac- 
ings of the fine lattice, a s . The fits (3.1C) 
are visualized by the dashed lines. 



where Si are two-point operators of the lattice mo- 
nopole charges corresponding to different separa- 
tions between the charges. The term Si corresponds 
to the zero distance between the lattice monopoles, 
$2 corresponds to the unit distance and so on (see 
ref. 17] for further details). 

The two-point coupling constants, fi, of the lat- 
tice monopole action are shown in figures |and| function the distance between the lat- 
tice points. The numerical data corresponds to lowest, T = 1.6T C , and highest, T = 4.8T C , 
available temperatures. The number of blocking steps is fixed to n s = 6 while the spatial 
spacings of the fine lattice, a s , are available for a set of values ranging from a s = .16<t -1 / 2 
till a s = .25a- 1 / 2 . 

the leading term in the 



According to eq. (2.38 
lattice monopole action for large lattice monopoles 
(b ^> Xd) must be proportional to the Coulomb 
interaction, 
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Figure 2: The same as in figure [l] but 
for the temperature T — 4.8T C . 



To check this prediction we fit the coupling constan- 
ts fi by the Coulomb interaction ( 3.10| ) treating Cc 
is the fitting parameter. The fits are visualized by 
the dashed lines in figures [j] and 0. 

As one can see from the figures, this one-para- 
metric fit works very good. The x 2 /d-o.f. parame- 
ter is of the order of unity for most fits while it is 

of order of 2 for the smallest lattice spacing, a^/a = 0.16. Note that the dependence of the 
lattice Coulomb interaction on the distance between the interaction points is not a mono- 
tonic function, as one can see from the fitting curves. Similar non-monotonic behaviour is 
also observed in the numerical data. 



4 
3 
2 
I 



■ T=1.92T, 
♦ T=1.4T, 



0.8 



0.6 



0.4 



0.8 



1.2 



1.4 



0.2 



5 TIT 



Figure 3: The pre-Coulomb coupling Cc 
and the fits of Cc by eq. ( 3.11 ) for various 
temperatures T. 



Figure 4: The product of the screening 
length and the continuum monopolc density, 
eq. (3^8), calculated from the lattice mono- 



pole action (in units of the string tension) . 



Using the fitting of the action we obtain the values of Cc for various sizes of the lattice 
monopole, by/a = .96... 1.5 at different temperatures, T = (1.6 . . . 4.8)T C . According to 
eq. (|2.38 ) the pre-Coulomb coefficient Cc(b,T) at sufficiently large lattice monopole sizes, 
b (b 3> Xd) should be as follows: 



C c (b,T) 



R(T) b 2 



(3.11) 



where R is the product of the screening length and the continuum monopole density, 
eq. (H). 

We present the data for the pre-Coulomb coefficient, Cc{b,T) and the corresponding 
one-parameter fits ( 3.11] ) in figure ^. The fit is one-parametric with R being the fitting 
parameter. Again we observe that the agreement between the data for Cc and the fits is 
very good. We show the quantity R vs. temperature in figure § 

The lattice Coulomb form of the action and proportionality of the pre-Coulomb term 
to b~ 2 at large values of b was established [^] also in the three-dimensional Georgi-Glashow 
model for the 't Hooft-Polyakov monopoles. These facts does not come unexpected from 
the point of view of the discussion above. 
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3.3 Lattice monopole density 

As we have seen from previous section the density of squared lattice monopole charge should 
also contain the information about the parameters of the Coulomb gas model. According 
to eq. (2.51) the large-6 asymptotics of the quantity (k 2 (b))/b 2 can be used to extract the 
product of the screening length and the continuum monopole density R, eq. ( |3.8p . We 
have measured the density of squared lattice monopole charge for all available tempera- 
tures and lattice monopole sizes. As an example we plot in figures [5| and ^ the quantity 
(k 2 (b))/b 2 vs. the lattice monopole size, b, for lowest and highest available temperatures, 
respectively. 
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Figure 5: The ratio (k 2 (b))/b 2 vs. lattice 
monopole size, b, at T = 1.6 T c . 



Figure 6: The same as in figure || but for 
T = 4.8 T c . 



Our theoretical expectations ( 2.51| ) indicate that the function (k 2 (b))/b 2 must vanish 
at small monopole sizes and tend to constant at large b. This behaviour can be observed in 
our numerical data, figure ||, up to some jumps for densities with different n s . We ascribe 
these jumps to the lattice artifacts emerged due to finiteness of the fine lattice spacing, a, 
and finite volume effects. 

According to eq. ( |2.51| ) we should know the 
large-6 asymptotics of (k 2 (b))/b 2 to get the quan- 
tity R, eq. (|3lj|). These asymptotics are approxi- 
mated by averaging of the appropriate n s = 8 data 
for which the behaviour of the function in ques- 
tion is almost flat. Then we get the product of 
the screening length and the continuum monopole 
density, R, depicted in figure 0. 

The result is very similar to the one obtained 
from the behavior of the lattice monopole action. 
However, the quantity R obtained from the lat- 
tice monopole density is a bit larger than the same 
quantity calculated from the lattice monopole ac- 
tion. This fact can be expected since we have ap- 
proximated the asymptotics of (k 2 (b))/b 2 by the average of the n s 
slightly differ from a correct asymptotics. 
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Figure 7: The same quantity as in fig- 
ure now calculated from the lattice 
monopole density. 
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3.4 Check of the Coulomb gas picture for continuum monopoles 

Although both quantities R p and R ac t correspond to the product of the screening length and 
the continuum monopole density, from a numerical point of view both i?'s are independent. 
To check the self-consistency of our approach we plot the ratio of these quantities in figure || 
It is clearly seen that the ratio is independent of the temperature and very close to unity, 
as expected. The 10%— 15% deviation of this quantity from unity may be explained by 
reasons mentioned in the end of the previous subsection. 
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Figure 8: Check of the self-consistency: the 
ratio of the quantities R, eq. ( |3.8| ) obtained 
from the lattice monopole action and density. 



Figure 9: Check of the dilute Coulomb gas 
picture: quantities C sp , eq. ( |3.12 ), calculated 
from the action and density. 



A check of the validity of the Coulomb gas picture can be obtained with the help of 
the quantity 



C sp (T) 



c. 



sp 



(T) 



<r sp (T) 



X D (T)p(T) R(T) 



(3.12) 



where a sp is the spatial string tension. In the abelian projection approach the spatial string 
tension should be saturated by the contributions from the static continuum monopoles. In 
the dilute Coulomb gas of monopoles the string tension is |p!8| : a = 8^/p/gM while the 
screening length is given by (^1]). These relations imply that in the dilute Coulomb gas of 
continuum monopoles we should get C sp = 8. 



We use the results for the spatial string tension of ref. [22] in the high temperature 
SU(2) gluodynamics. It was found that for the temperatures higher than T ~ 2T C the 
spatial string tension can be well described by the formula: a sp {T) = 0.136(11) g\ D (T) T 2 , 
where <74_d(T) is the four-dimensional SU(2) 2-loop running coupling constant, 



11 



12vr 



2 l0g (A^)+44^2 l0 S 



2l0g (B 



with the scale parameter At = 0.076(13) T c . Taking also into account the relation between 
the critical temperature and the zero-temperature string tension |^], T c = 0.69(2) y/a we 
calculate the quantity C sp and plot it in figure |9| as a function of the temperature, T. If 
the Coulomb picture works then C sp should be close to 8. Prom figure ^ we conclude that 
this is indeed the case at sufficiently high temperatures, T/T c > 2.5. 
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4. Conclusions 

In order to describe the lattice monopole dynamics we have proposed to consider the lattice 
monopoles as the defects which are blocked from continuum. In other words the lattice 
was suggested to be a measuring device for the continuum monopoles. As a result we are 
able to draw the following conclusions: 
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Using the Monte Carlo results for the density of the squared monopole charges and 
the monopole action we are able to calculate the product of the abelian magnetic 
screening length and the monopole density corresponding to the continuum Coulomb 
gas model. The values of this parameter obtained from density and action measure- 
ments are consistent with each other and — at sufficiently high temperatures — are 
consistent with known results for the spatial string tension. At temperatures T j 
T c > 2.5 the spatial string tension is dominated by contributions from the continuum 
static monopoles. 

The continuum Coulomb gas model can describe the results of the Monte Carlo 
simulations for the action and density of the lattice monopoles. The dependence of 
these quantities on the physical sizes of the lattice monopoles (the size of the cell 
which is used for the monopole detection) is in a good agreement with the analytical 
predictions. 

The lattice monopole action is dominated by the mass and the Coulomb terms for, 
respectively, small and large sizes of the lattice monopoles. A relation between the 
density of the squared magnetic charges and the monopole action is established 
(eqs. ( fTl6D , ( ET47D , flggp and/or eqs. ( pj >, (|236|) , (gjg)). Our analytical deriva- 
tion was done in a gaussian approximation. We have shown that the corrections to 
the leading terms of the action are small in the large-& limit. A detailed analysis of 
the corrections will be discussed elsewhere |P 



We believe that this method can also be applied to the four dimensional non-abelian 
gauge theory. We think that this would allow to get (at least, a part of) parameters of the 
dual superconductor model. 
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